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Abstract 

Our theoretical comparisons suggest that the ATF measurement 
results of April 2000 for energy spread, bunch length, and hori- 
zontal emittance vs. current, and a low current emittance ratio 
of about 1% are generally consistent with intrabeam scattering 
(IBS) theory, though the measured effects appear to be stronger 
than theory. In particular, the measured factor of 3 growth in 
vertical emittance at 3 mA does not seem to be supported. It 
appears that either (1) there is another, unknown force in ad- 
dition to IBS causing emittance growth in the ATF; or (2) the 
factor of 3 vertical emittance growth is not real, and our other 
discrepancies are due to the approximate nature of IBS theory. 
Our results suggest, in addition, that, even if IBS theory has 
inaccuracies, the effect will be useful as a diagnostic in the ATF. 
For example, by measuring the energy spread, one will be able to 
obtain the emittances. Before this can be done, however, more 
benchmarking measurements will be needed. 
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Intrabeam Scattering Analysis of ATF Beam Data 
Taken in April 2000 

K.L.F. Bane, H. Hayano, K. Kubo, T. Naito, T. Okugi, and J. Urakawa 
Introduction 

Between April 13-19, 2000 single bunch energy spread, bunch length, and 
horizontal and vertical emittances in the ATF were all measured as a func- 
tion of current jjj- 1| . One surprising outcome was that the vertical emittance 
appeared to grow by a factor of 3 by a current of 3 mA. The ATF is a proto- 
type damping ring for the JLC/NLC linear colliders, and the concern with 
this result is that it may portend an as yet not understood and unexpected 
growth in such damping rings, which would have negative ramifications on 
collider performance. However, since the x-y coupling in the ATF is very 
small (~ 1%), and since the emittance measurements were performed on 
the beam after it had been extracted from the ring, the question was, How 
much of this measured y-emittance growth was real and how much was due 
to measurement error, such as dispersion in the extraction line or in the wire 
monitors used for the measurements. 

With the ATF as it is now, running below design energy and with the 
wigglers turned off, the beam properties are strongly affected by intra-beam 
scattering (IBS), an effect that couples the three dimensions of the beam 
together. In April 2000 all the beam dimensions were measured to varying 
degrees of accuracy, and the hope is that the knowledge of IBS theory can 
be used to check for consistency in the data. Besides the question of the 
vertical emittance growth, we hope that IBS theory can be used in beam 
diagnostics in the future. For example, the energy spread measurement 
is quick, accurate, and easy to do. It would be nice if this measurement 
could be used to estimate the beam emittance and/or bunch length, beam 
properties that are much more difficult to measure directly. 

The literature on intrabeam scattering is quite extensive. (For an in- 
troduction, see for example the IBS section and its bibliography, written 
by A. Piwinski, in the Handbook of Accelerator Physics and Engineering 
[Q]). The first rather thorough treatment of IBS in circular accelerators is 
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due to Piwinski (P), derived following a two-particle Coulomb scattering 
formalism^]. Another formalism is that of Bjorken-Mtingwa (B-M), ob- 
tained following quantum mechanical two particle scattering rules ||. This 
is the formalism that is more often used in modern optics programs that also 
calculate IBS, programs such as SAD[f7|. The B-M result is considered to be 
more general in that the combination of optics terms /3 x (j) x = {(3 x rj' x — P' x rj x /2) 
around the ring do not need to be small compared to rj x , whereas in the P 
method it seems they do|J. (Note that this condition is typically violated 
in modern low emittance storage rings.) The B-M formalism, however, does 
not include vertical dispersion, whereas the P formalism does. Neither for- 
malism includes x-y coupling, though a more generalized formulation, which 
includes both linear coupling and can also be applied to low emittance ma- 
chines, is given by Piwinski in Ref. ||. Note that in deriving such for- 
mulations many approximations were made, having to do with the cut-off 
distance for scattering (typically taken to be the vertical beam size), curved 
trajectories, etc. In addition, it appears that no current IBS formalism prop- 
erly accounts for the effects of potential well distortion or of the micro-wave 
instability. Finally, T. Raubenheimer has pointed out that IBS does not re- 
sult in Gaussian bunch distributions, and that these theories give rms beam 
sizes that can greatly overemphasize few particles occupying the tails 1 1C]. 

Most of the early papers comparing IBS theory and measurement were 
for bunched and unbunched hadronic machines, where the effect tends to 
be more pronounced than in high energy electron machines. For example, 
Conte and Martini, for unbunched protrons in the CERN Antiprotron Ac- 
cumulator ring, found good agreement for the longitudinal and horizontal 
IBS growth rates, but no agreement for the vertical rate|ll|. Evans and 
Gareyte, for bunched protrons in the SPS, found good agreement for the 
radial emittance growth rate with time, once an additional factor repre- 
senting gas scattering was included[12|. As for electron machines, both IBS 
and the related Touschek effect have become important effects in modern, 
low emittance light sources. C.H. Kim at the ALS found that he can get 
agreement with measured horizontal emittance growth with current, but to 
do this he needed to include a significant additional fitting factor in the 
calculations fi"3|1 . So it may be too much to expect good agreement between 
IBS theory and measurement without the use of such fudge factors. Indeed, 
A. Piwinski has said that, given the approximations taken in deriving IBS 
theory, one can expect agreement between theory and measurement only on 
the order of a factor of 2 [|| . Yet even if this became the case, it may still be 
possible to benchmark the ATF with accurate emittance and energy spread 
measurements, to find the fudge factors. Then in the future, one may be able 
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to perform simpler measurements — like the energy spread measurement — to 
get an estimate of parameters that are more difficult to obtain directly, such 
as bunch length and emittance. 



Piwinski's Solution 

We are interested in what happens at low coupling, and we will concentrate 
on using the more simplified version of Piwinski's solution. We begin by re- 
producing the Piwinski solution in its entirety 0. Note that there is nothing 
new in this section, except the way potential well distortion is added to the 
calculation. 

Let us consider the IBS growth rates in energy p, in the horizontal di- 
rection x, and in the vertical direction y to be defined as 
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Here a p is the rms (relative) energy spread, e x the horizontal emittance, and 
e y the vertical emittance. According to Piwinski the IBS growth rates are 
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where the brackets () mean that the enclosed quantities, combinations of 
beam parameters and lattice properties, are averaged around the entire ring. 
Parameters are: 
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/ (a,M)=8^ 1 {21„[|(I + ; L)]-0. 5 77...} 



du (8) 



The function / is given by: 

l-3n 2 
PQ 

p2 = a 2 + (l-a 2 )u 2 , Q 2 = b 2 + (1 - b 2 )u 2 (9) 

The global beam properties that are affected by IBS are the rms (relative) 
energy spread a p and the rms bunch length a s , the horizontal emittance 
e x , and the vertical emittance e y . Other global properties are the bunch 
population N, the relative velocity /3, and the energy factor 7. Note that ro 
is the classical radius of the particles (for electrons ro = 2.82 x 10 -15 m). The 
lattice functions needed are the beta functions (3 X , @ y , and the dispersion 
functions D x , D y . Note also that a x p = y/f3 x e x , and a y p = y/(3 y e y . The 
parameter d represents a cut-off for the IBS force, which Piwinski says should 
be taken as the vertical beam size, but he also points out that the results 
are not very sensitive to exactly what is chosen for this parameter. 

Then the steady-state beam properties are given by 

€x = f£2 , €y = € Jfi ; CT 2 = *po ^ (1Q) 

1 T x /T x 1 T y /T y 1 Tp/Tp 

where subscript represents the beam property due to synchrotron radia- 
tion alone, i.e. in the absence of IBS, and t x , T y , and t p are the synchrotron 
radiation damping times in the three directions. These are 3 coupled equa- 
tions in that all 3 IBS rise times depend on 63;, 6y, and dp. Note that a 4th 
equation, the relation between a s and a p , is also implied; generally this is 
taken to be the nominal (zero current) relationship. 

Piwinski suggests iterating Eqs. until a self-consistent solution is 
found. Our experience is that this has the problem that negative values 
of emittance or a 2 can be obtained with this procedure, causing difficulty 
in knowing how to continue the iteration. We find that a better method is 
to convert these equations into 3 coupled differential equations, such as is 
done in Ref. |[g| ] to obtain the time development of the beam properties in 
a ring. Here, however, we use it only as a mathematical device for finding 
the steady-state solution. For example, the equation for e x becomes 

de x _ _{e x — e x o) ^ e x 

and there are corresponding equations for e y and a 2 . 

We have three comments: (1) In a storage ring vertical emittance is 
generally the result of two phenomena, x-y coupling — caused, for example, 
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by rolled magnets — and vertical dispersion — caused by vertical closed orbit 
distortion. Although our formalism technically is valid only for the case of 
no x-y coupling, we believe that it can also be used for the case of weak 
coupling. In such a case we simply pick a zero-current emittance ratio 
r xyO = fj/oAzO) an d then set e y o accordingly. Note that r xy0 is meant to 
include both the contributions of x-y coupling and of vertical dispersion. 
Finally, note that the vertical emittance, assuming only vertical dispersion, 



can be approximated by [14] 



'y) rms _z 



with J e the energy damping partition number. 

(2) It is not clear how impedance effects and IBS effects interact. At the 
ATF we have shown that up to the currently highest attainable single bunch 
currents (~ 3 mA) a micro- wave threshold has not yet been reached Q. But 
we do appear to have a sizeable potential well bunch lengthening effect. 
We will assume that this case can be approximated by adding the proper 
multiplicative factor f pw {N) } obtained from measurements, to the equation 
giving a s in terms of a p . (3) As mentioned earlier, the IBS bunch distri- 
butions are not Gaussian, and tail particles can be overemphasized in the 
beam size solutions given above. T. Raubenheimer in Ref. |n| gives a way 
of estimating IBS beam sizes that better reflect the particles in the core of 
the beam. Note that the optics computer program SAD can solve the IBS 
equations of B-M, include the effects of orbit errors, and also find the sizes 
of the core of the beam. 



Parameter Studies 

We have programmed the above Piwinski equations. We have also pro- 
grammed the B-M equations, which also give 3 IBS growth rates (though 
there is a factor of 2 difference in their definition). Let us begin our nu- 
merical studies by comparing the results of the two programs when applied 
to the ATF lattice and beam properties. (Note that there is a published 
comparison of results of these two methods applied to the CERN AA ring 



in Ref. IE]; however, the variation of the lattice parameters around the ring 
was not included in the P method calculation.) As parameters we take: 
E = 1.28 GeV, a p0 = 5.44 x 10" 4 , a s0 = 5.06 mm (for an rf voltage of 
300 kV), e X Q = 1.05 nm, t p = 20.9 ms, t x = 18.2 ms, and t x = 29.2 ms. 
The ATF circumference is 138 m, {(3 X ) = 4.2 m, (/3 y ) = 4.6 m, (D x ) = 5 cm. 
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The function /3 Z |(/> X | is about .5D X at positions of the minima in D x , and 
about .15D X at positions of the maxima, which are not so small, and we 
might expect some disagreement between the two methods. Note that for 
the averages represented by brackets in Eqs. [2||4|, and their counterparts in 
the B-M method, we calculate, as is normally done, the appropriate com- 
bination of lattice and beam properties at the ends of the lattice elements, 
connect these with straight lines, and then find the average of the resulting 
curve. 

Fig. H displays the 3 differential IBS growth rates as obtained by the 
two methods (blue for P, red for B-M) when / = 3.1 mA and r xy Q = .01. 
The IBS growth rates are just the average values over the ring of these 
functions. Since the ATF has two- fold symmetry we give the result for only 
half the ring (the straight section is in the middle). Here we include no 
vertical dispersion, i.e. the vertical emittance is presumed to be only due to 
(weak) x-y coupling, and f pw is set to 1. We see almost perfect agreement 
for the differential growth rates obtained by the two methods, with only 
slight differences in the peaks. As for the averages, for the Piwinski method 
l/T p = 28.4 s"\ 1/T X = 24.8 s _1 , and 1/T y = .67 s -1 . The B-M results 
are 28.9, 26.2, and .58 s _1 , respectively. Note that the y growth rates 
are very small. As for the steady state beam properties, for this example 
Cp/cpO = 1-566, e x /e x o = 1.825, and e y /e x o = .0102 for P, and 1.587, 1.910, 
.0102 for B-M, respectively. We see that for the ATF the results of the 
two methods are almost the same. Finally, note that when making the same 
comparison for the ALS at Berkeley, a low emittance light source with /3 x \<f>x\ 
comparable to D x , we again find very close agreement in the results of the 
two methods. 

In the ATF the rms vertical dispersion, after correction, is typically 3- 
4 mm. To simulate this effect we added a randomly generated component 
of D y (weighted by y/~j3y~) at the high and low @ y points, and connected 
these with straight lines (see Fig. ^a). This was added to the P calculation. 
In the ATF typically r xy0 » .01 and J e = 1.4. We see from Eq. |TJ that if 
(Dy) r ms = 4 mm, then 30% of the low current emittance is due to dispersion, 
and the rest due to x-y coupling. Let us consider an example now where 
the size of e y o is given entirely by vertical dispersion, with r xy o = .01 and 
(D y )rms = 7.5 mm. In Fig. ||b we plot the resulting 5{1/T y ) at I = 3.1 mA. 
The growth is much larger than before, with 1/T y = 11.8 s _1 . (Note that 
Cp/(TpQ = 1.48, e x /e x o = 1.71, and e y /e x o = .015.) We can understand the 
change in 1/T y if we look back at the equations for the growth rates in x 
and y, Eqs. |3| and ||. In Eq. ^ the horizontal growth rate is dominated by 
the second term in the brackets, the one proportional to D^/a^. In Eq. || 
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Figure 1: Differential growth rates over 1/2 the ATF, as obtained by the 
Piwinski (blue) and the Bjorken-Mtingwa (red) methods Here I = 3.1 mA, 
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Figure 2: Random vertical dispersion [(D y ) rms = 7.5 mm] (a), and the 
resulting differential vertical growth rate (b) over 1/2 the ATF. X/T y = 
11.8 s _1 . Here I = 3.1 mA and r xy o = .01. 
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the vertical growth rate, when D y = 0, is given by the small first term in 
the brackets. Since the second term is proportional to Dy/oZp, 1/T y will 
become comparable to 1/T X when (D y ) rrns ~ {D x )^/e y /e x (since (f3 x ) and 
{(5 y ) are similar), which equals 5 mm. Note that the importance grows as 
the second power of (D y ) rms . 

We see that the Piwinski and the Bjorken-Mtingwa methods give essen- 
tially the same solution for the ATF, and since the P method allows for 
vertical dispersion, we will choose to continue our simulations using this 
method. In Figs. ||||5] we give the steady-state emittance and energy spread 
vs. various parameters in the ATF according to Piwinski's IBS theory. 
We include curves representing (D y ) rms = 0, 3, and 6 mm. [Note that 
{D y ) r ms = is not a realistic condition for the ATF.] Given the approxi- 
mate nature of IBS theory, these results are meant to give an idea of the 
sensitivities of the steady-state beam properties to various parameters, and 
not to give absolutely correct predictions. In Fig. || we show the depen- 
dence on I, with f pw = 1, r xy o = .01. We note that the vertical emittance 
growth is almost zero with zero vertical dispersion. At {D y ) rms = 6 mm, 
e y has increased by 44% by I = 3 mA. In Fig. |4] we show the dependence 
on r xy o, with f pw = 1, / = 3.1 mA. We note that at r xy o = the vertical 
emittance does not go to zero when (D y ) rms is not zero. Note also that for 
the case (D y ) rms = 6 mm, the a p (r xy o) curve is rather linear over the range 
shown, with a slope Aa p /a p / Ar xy o = 6.6; i.e. an (absolute) change in r xy o 
of 0.005 produces a relative change in a p of only 3%. In Fig. § we give the 
dependence on f pw , with I = 3.1 mA, r xy Q = .01. Note that at I = 3 mA, 
the ATF measurements give f pw = 1.25. We see that the energy spread and 
emittances are not very sensitive to this parameter. For example, a change 
of .1 in f pw yields a 1% change in a p . 

Measurements [[Ij] - 13|] 

In April 2000, a p , a s , e x , and e y were all measured as functions of current to 
varying degrees of accuracy. Among the most accurate is believed to be the 
energy spread measurement, which is performed on a thin screen at a high 
dispersion region in the extraction line (D x = 1.73 m). At different currents 
the measured beam width was fit to a Gaussian (the fits were very good) and 
the rms width was extracted. This measurement was performed on April 14 
(see Fig. ^a). The rms scatter in the extracted a p was less than 2%, and the 
precision should be better than 1%. Note that although data was obtained 
for currents up to ~ 2 mA only, from experience we have confidence that we 
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Figure 3: Calculations of the current dependence of the beam properties 
when r xy0 = .01, f pw = 1. 
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Figure 4: Calculations of the r xy o dependence of the beam properties when 
/ = 3.1 mA, f pw = 1. 
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Figure 5: Calculations of the f pw dependence of the beam properties when 
I = 3.1 mA, r xyQ = .01. 




Figure 6: Measurements of energy spread (a) and bunch length (b), with 
V c = 300 kV, performed on April 13-14, 2000. The fits (the curves) are 
chosen to give the expected zero current results. 
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can extrapolate to ~ 3 niA. 

The rms bunch length, using a streak camera, was also measured about 
the same time, and there was more scatter in the data(see Fig. ^b). With 
streak camera measurements, however, there is always the question of whe- 
ther space charge in the streak camera itself could have added a systematic 
error to the results. Checks were made with light filters, so we don't think 
this is a problem, but this still adds a slight uncertainty to the results. Both 
results were fit to a smooth curve, chosen to give the expected zero current 
results (see Fig. |6|). Note that these results, if true, imply an extremely large 
bunch length increase at low currents. 

The emittances were measured on wire monitors in the extraction line on 
April 19. The results are reproduced by the plotting symbols in Fig. [7|b-c. 
We see that the x emittance appears to grow by ~ 80% by I = 3 mA; the y 
emittance begins at about 1% of the x emittance, and then grows by a factor 
of 3. We believe that the x-emittance measured should be fairly accurate. 
The y-emittance, however, since it is so small, could be corrupted by many 
factors, such as dispersion in the extraction line or the wire monitors (roll of 
the measurement wires, for example, has been checked and shown not to be 
significant). Note that it appears that r xy o ~ .01. We estimated above, using 
Eq. that if the vertical emittance is dominated by vertical dispersion, 
then r xy o = .01 implies that (D y ) rms = 7.5 mm, which is significantly larger 
than the measured 3-4 mm. Unfortunately, we do not know what D y (s) was 
during the April measurements, and therefore cannot make a more precise 
comparison with calculation. 

Note that all the measurements were not performed on the same day, 
and since IBS depends on the status of the machine (e.g. on the vertical 
closed orbit distortion), it is possible that the longitudinal and transverse 
measurements correspond to slightly different machines as far as IBS is con- 
cerned. For example, on April 13, the day before the a p measurements of 
Fig. H the energy spread was also measured. Those results, when fitted and 
extrapolated to / = 3 mA, gave a 2% smaller rms value than given here. 
This adds some uncertainty to our results. 

Comparison with IBS Theory 

To study the consistency of these measurements with IBS theory, we per- 
form IBS calculations where we take the potential well factor f pw that was 
measured, and adjust r xy o to obtain the measured a p (and a s ) at / = 3 mA, 
for the cases (D y ) rms = 3 and 6 mm. [Remember: the typical measured 
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value is (D y ) rms = 3-4 mm.] The best fits were found for r xy o = .0104 and 
.006, for the cases (D y ) rms = 3 and 6 mm, respectively. The results are 
shown in Fig. |^, where they are compared to the measured data. (Note that 
in Fig. 0a, we reproduce, using plotting symbols, the smooth curve fits to 
the measured data of Fig. ||.) 

We first notice that whether there is a 3 mm or 6 mm y dispersion the 
fitted results are very similar. We find for our fitted results that the e x de- 
pendence on current is in reasonable agreement with measurement (though 
the calculated result is low by 10% at 3 mA), and that r xy o ~ .006-. 01, 
which is a little low compared to .01. The biggest discrepancy, however, 
seems to be that the vertical emittance growth is much lower for the calcu- 
lations than the measurements: for the case (D y ) rms = 3 mm we see 10% 
growth by 3 mA; for {D y ) rms = 6 mm we see 60% growth by 3 mA, which, 
however, is not close to the factor of 3 that was measured. 
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Figure 7: ATF measurement data (symbols) and IBS theory fits (the 
curves). The symbols in (a) give the smooth curve fits to the measured 
data of Fig. ||. In the theory the parameter r xy o was adjusted to get the 
measured a p at / = 3 mA. 

Let us suppose that the measured vertical emittance growth is not real, 
and is due to measurement error. In such a case, our results appear to give 
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reasonable agreement between IBS theory and the ATF measurements, for 
a p , a s , and e x vs. I, and low current e y , all without extra fudge factors. For 
the present set of data, however, we are still left with some uncertainties. 
For example, it was suggested earlier that the accurate knowledge of a p {I) 
is important. It was shown that for a realistic type of vertical dispersion 
for the ATF, {D y ) rms = 6 mm, at / = 3.1 mA, a change in r xy Q of .005 
produces a relative change in o p of only 3%. Or, conversely, if the inaccuracy 
in <7 P (3 mA) is 3%, then the fitted r xy Q shifts by .005, which is not small 
compared to .01. And from the measured difference in a p on April 13 and 
14, it seems that the uncertainty in a p is of this order. Finally, we should 
note that another source of uncertainty is that we don't know the function 
D y (s) on the date of the measurements. 

However, there is one major quantative discrepancy between theory and 
measurement, which has to do with non-Gaussian beam tails. SAD finds 
that for the ATF the IBS induced emittance growth, when not counting such 
tail particles, is only 2/3 of that when all particles are included (or 80% in 
the case of a p ), in agreement with calculations for the ATF given in Ref. [1C]. 
Our calculations include all beam particles, while the measurements consider 
only core particles. If the tail particles are indeed as significant as these 
results suggest, then the effect measured at the ATF is much stronger than 
predicted by IBS theory, and there would no longer be good agreement with 
measurement. If one wants to think about a fudge factor that is multiplied 
with the bunch charge, then we estimate that this factor must be ~ 1.5-2.0 
(though no single such factor gives good agreement for all cr p (I), cr s (I), and 
e x (I)). From this result, it appears that either there is another, unknown 
force in addition to IBS causing emittance growth in the ATF, or the problem 
is the approximate nature of IBS theory. 

In spite of this discrepancy, we believe that the IBS effect will be useful 
as a diagnostic at the ATF. However, before we can use it as such, we need 
to perform more benchmarking experiments. Such measurements include 
the effects of varying the rf voltage, the vertical dispersion, and the x-y 
coupling. Additional independent measurements are also desireable, such 
as using the interferometer to measure the beam sizes in the ring. The 
calculations can also be improved by including a cut-off for the tails. In 
addition, one might, for example, include the measured vertical dispersion 
(instead of a randomly generated one) and/or use Piwinski's more involved 
formulation that properly includes linear coupling. 
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Conclusion 



Our theoretical comparisons suggest that the ATF measurement results of 
April 2000 for energy spread, bunch length, and horizontal emittance vs. 
current, and a low current emittance ratio of about 1% are generally con- 
sistent with intrabeam scattering (IBS) theory, though the measured effects 
appear to be stronger than theory. In particular, the measured factor of 
3 growth in vertical emittance at 3 mA does not seem to be supported. 
It appears that either (1) there is another, unknown force in addition to 
IBS causing emittance growth in the ATF; or (2) the factor of 3 vertical 
emittance growth is not real, and our other discrepancies are due to the ap- 
proximate nature of IBS theory. Our results suggest, in addition, that, even 
if IBS theory has inaccuracies, the effect will be useful as a diagnostic in the 
ATF. For example, by measuring the energy spread, one will be able to ob- 
tain the emittances. Before this can be done, however, more benchmarking 
measurements will be needed. 
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